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^ ■ 1 Introduction 
cn 

\ Solving both classical and quantum Hamiltonian systems follows the pattern: the more "struc- 

■ ture" (i.e., symmetries, first integrals, etc.) a Hamiltonian system has the more likely it can be 

CN ■ solved explicitly. Superintegrable systems are very rich in structure, which makes it possible to 

solve them explicitly both analytically and algebraically. In this paper we investigate maximally- 
superintegrable Hamiltonian systems of classical mechanics defined by natural Hamiltonians 
^ . which admit quadratic in the momenta first integrals of motion. Recall that a Hamiltonian 

?H \ system with n degrees of freedom is said to be maximally-superintegrable if it admits 2n — 1 

global, functionally independent, single-valued integrals of motion that do not depend explicitly 
on time. In what follows we also assume that such systems admit quadratic in the momenta 
first integrals of motion, which under some additional conditions afford orthogonal separation 
of variables in the associated Hamilton- Jacobi equation. In fact, such superintegrable systems 
can be solved via orthogonal separation of variables in more than one way, which is the reason 
they are called multi-separable. This particular class of superintegrable systems includes such 
important examples as the Kepler potential, Calogero-Moser model, harmonic oscillator and 
others. 

Many efforts have been made over the years to classify superintegrable Hamiltonian systems 
defined in spaces of constant curvature of the type described above. Formally, it was initiated 
in 1965 by Fris et al who produced [10] a list of superintegrable potentials of this type defined 
in Euclidean plane. In 1990 Evans [7] extended this classification to superintegrable potentials 
defined in Euclidean 3D-space. More recently, several authors have launched the project of 
classification of superintegrable systems, based on the algebraic properties of the symmetry 
algebra associated with a given superintegrable system. This idea has proven to be quite fruitful 
since the symmetry algebra of a superintegrable Hamiltonian system has a rich structure, unlike 
the symmetry algebra of a completely integrable Hamiltonian system admitting only n first 
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integrals in involution, whose symmetry algebra is merely abelian (for more details refer to 
Kalnins et al |16j and the relevant references therein). 

Another approach to the problem of classification of superintegrable systems of the type 
described above defined in the Euclidean plane has been realized by Adlam et al [5] within 
the framework of the invariant theory of Killing tensors (ITKT) (see [H [21 [3l IH O [H [TTl [121 
[13 [13 [m [M [201 [221 [251 [221 [27j and the relevant references therein for more details). More 
specifically, the authors derived an invariant characterization of the superintegrable potential of 
the Kepler problem, using invariants and joint-invariants of Killing two-tensors that determined 
first integrals admitted by the system in question. 

In this work we extend some of the ideas employed in [2] to characterize, using joint invari- 
ants of Killing two-tensors [22] defined in the Euclidean plane, the well-known Smorodinsky- 
Winternitz (SW) potential and explain the geometric meaning of its arbitary parameters. We 
also verify that the important integrable perturbation of the Smorodinsky-Winternitz poten- 
tial known as Tremblay-Turbiner-Winternitz (TTW) potential loses the property of multi- 
separability for all but one value of the arbitrary parameter k (refer to the formula ([7])). 

2 Superintegrable, multi-separable systems in Euclidean plane 

Let (A^,g) be an n-dimensional pseudo-Riemannian space of constant curvature. Consider a 
Hamiltonian system defined in the cotangent bundle T*M by the following natural Hamiltonian 
function 



H{q,p) = ^g'HQ)PiPj + y{Q), hj = h---^n, (1) 

where q = {q^, . . . , g"), p = (pi, . . . ,Pn) are the corresponding canonical coordinates, g^^ denotes 
the components of the metric tensor, and V{q) is the potential function. In what follows we will 
study the Hamiltonian systems defined by ([T]) that admit 2n — 1 first integrals of motion, given 
by single- valued, functionally independent, globally defined functions 

Fi = /f, F2, . . . , F„, F„+i, . . . , F2n-i G T{T*M). (2) 

that enjoy the property 

{Fi,Fj} = {Fk,Fe} = 0, i,j = l,...,n;k,i = l,n + l,...,2n-l, (3) 

where {•, •} denotes the corresponding (canonical) Poisson bracket defined in J^{T*Ai). A Hamil- 
tonian system described above is said to be a maximally- superintegrable Hamiltonian system 
defined by a natural Hamiltonian ([T|). Furthermore, let the first integrals Q be quadratic in 
the momenta p that is of the form 

F,{q,p)=Ki'{q)p,p,+U,{q), (4) 

where i = 1, . . . ,2n — 1, i, j = 1, . . . ,n. As is well known, the condition {H, Fi} = implies for 
each £ = 1, . . . ,2n — 1 that 

[Ke,g]=0, (5) 

where [•, •] denotes the Schouten bracket [2l] and 



dikedV) = 0, 



(6) 
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where Ki = Kig~^. Equation ([5]) is called the Killing tensor equation, implying that K^'^ are 
the components of a Killing two-tensor Ki, while Equation ^ is said to be the correspodning 
compatibility condition (known also as the Bertrand-Darboux equation(s)). 

Naturally, one of the most important problems that concerns this type of superintegrable 
Hamiltonian systems is the problem of classification. Thus far it has been shown that it can be 
approached from different directions by taking into account the algebraic, analytic and geometric 
properties of superintegrable potentials in question and their corresponding first integrals of 
motion. 

First of all, different superintegrable potentials can be derived once the Killing tensors defining 
the quadratic part of first integrals of motion Q have been classified (see Section [2D by solving 
the compatibility condition ([6]) for various sets of Killing two-tensors that determine quadratic 
in the momenta first integrals of motions (jl]). This approach has been successfully employed, 
for example, to derive lists of superintegrable potentials defined in the Euclidean plane and 
3D-space (see [10] and [7] respectively) 

Next, as it has been shown by various authors in recent years (see [16] and the relevant 
references therein) one can use in this context the structrue algebra of first integrals associated 
with a given superintegrable Hamiltonian system in view of the fact that it is not abelian ([3|). 
More specifically, for a given superintegrable Hamiltonian potential ([1]) that admits 2n — 1 first 
integrals of motion ([2|) one can determine the order at which the algebra of first integrals closes 
with respect to the Poisson bracket, which is an algebraic characteristic of the superintegrable 
potential in question. 

Alternatively, one can make use of the invariant properties of the Killing two-tensors that 
appear in ^ by computing their joint invariants [22] and characterizing the associated super- 
integrable potential accordingly. This approach has been employed by Adlam et al [2] to derive 
an invariant characterization of the Kepler potential defined in the Euclidean plane ¥?. 

In this work our goal is twofold. First, we wish to characterize the Smorodinsky-Winternitz 
(SW) superintegrable system defined in in terms of Cartesian coordintes {q^^g^) = {x,y) by 
the following natural Hamiltonian function: 



H = \{pl+pl)-u:{x' + y') + ^,+^,. (7) 

This system is an integrable perturbation of the harmonic oscillator; its potential was derived 
in [To] in the framework of the problem of classification of superintegrable potentials defined in 
¥? that admitted two quadratic in the momenta first integrals of motions ([3]), using the first 
approach reviewed above. In what follows we target to derive an invariant characterization of 
the potential of ([7]) in terms of the joint invariants of the associated Killing two tensors and 
explain the meaning of the arbitrary parameters w, a and (3 from this viewpoint. 

Our second goal is to investigate a generalization of the SW potential given by ([7]) and known 
as the Tremblay-Turbiner-Winternitz (TTW) potential. The TTW potential was obtained by 
Trembley et al [231 [24] by rewriting the quantum analogue of the natural Hamiltonian ([7|) in 
polar coordinates (r, 9) and then replacing 9 with kO, thus introducing a new parameter k. The 
resulting natural Hamiltonian 



H 
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cos^ k9 sin^ k9 



(8) 



is characterized by a potential that preserves some of the properties of the SW potential de- 
termined by ([7]). In particular, it is obviously still separable in polar coordinates. The big 
question is whether or not the system is still superintegrable. The authors conjectured that 
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the system defined by the potential of ([8]) was classically and quantum superintegrable for all 
rational values of k. This conjecture has had a major impact on the development of the theory 
of superintegrable systems, resulting in a host of researchers trying to prove it (for more details, 
see [16] and the relevant references therein). In this work we will show that the Hamiltonian 
system defined by ([8]) is superintegrable in the sense of admitting quadratic first integrals of 
motion iff k = ±1, that is when it is of the form defined by ([7]). 

3 Invariant theory of Killing tensors in Euclidean plane 

Recall that the invariant theory of Killing tensors (ITKT) has been developed to study intrinsic 
properties of Killing tensors defined in pseudo-Riemannian spaces of constant curvature modulo 
the action of the corresponding isometry groups. This arrangement makes it an analog of the 
classical theory of invariants of homogenous polynomials which are studied within the framework 
of the latter theory under the action of the general linear group or its subgroups (see Olver [18] 
and the relevant references therein for more details). 

The origins of the theory can be traced back to the 1965 paper by Winternitz and Fris [25] 
in which the authors studied second order symmetries of the Laplace operator defined in 
modulo the action of the special Euclidean group SE{2). This problem is the quantum analog 
of the corresponding problem in classical mechanics involving the Killing two-tensors defined 
in under the action of SE{2), which was independently solved in 2002 by McLenaghan et 
al [17]. In the latter case the authors not only computed the invariants of the group action 
5^(2) O K?{¥?), where K?{¥?) denotes the vector space of Killing two-tensors in E^ , and used 
them to discriminate between the different orbits of the group action, but also computed in each 
case the corresponding moving frames map [HUH [18], as well as applied the resulting algorithm 
to solve Hamiltonian systems defined on E^ that arised in physics via separation of variables. 
These results laid the groundwork that allowed a number of authors to successfully employ 
this approach in a more general setting to study Killing tensors defined in spaces of constant 
curvature (for more details, see [B [11 [3l ll [3 [111 [13 [111 [III [IS [201 [22l [sll [2^ 

relevant references therein). 

In what follows our notations are compatible with those adopted in [2]. In this section 
we adapt the ITKT to derive in Section [J] an invariant characterization of the Smorodinsky- 
Winternitz potential given by ([7|). Since our goal is to characterize a superintegrable system 
defined in E^ that admit two distinct quadratic first integrals (other than the Hamiltonian), we 
first discuss how the theory works when there are two Killing two-tensors that define the first 
integrals of the form ^ . Consider the product space P defined by 

P:=/C2(e2) x/C2(e2), (9) 

in which each of the two factors is given by the following formula: 

K = (/3i + 2/34y + /?6y')5i 9i 

+(/33 - Pax - /35y - /3exy)di 82 (10) 

+ (/32 + 2/35X + /36x')92 52, 

representing the solution space of the Killing tensor equation ([5|) for the case = E^, where 
the arbitrary six parameters Pi, i = 1, ... ,6 are the corresponding constants of integration and 
di = ^,d2 = ^. Note dim/C2(E2) = 6. 

Recall that the action SE{2) O E^ given by 



X = X cosp3 — y sinp3 -I- pi, 
y = X sinps -|- y cosps + p2 



(11) 
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induces the corresponding action SE{2) O P'- 



di 


= ai cos^ ps — 2a3 cos sin + a2 sin^ — 2p2 cos 




-2^205 sinp3 + aepl, 


0.2 


= ai sin^ p3 — 2a3 cosp3 sinp3 + 02 cos^ ps — 2pia5 cosp3 




+2pia4 sinp3 + aepf, 


"3 


= (ai - a2)sinp3Cosp3 + a3(cos^p3 - sin^p3) 




+(pia4 +P2a5)cosp3 + (pias -p2a4)sinp3 - Q6P1P2, 


a4 


= Q4COSP3 + a5sinp3 - a6P2, 


"5 


= 05 cos P3 — 04 sin p3 — oqPi, 




= ae, 




= /^i cos^p3 - 2/^3 cosp3 sinp3 + /:i2 sm P3 - 2p2P4 cosp3 




-2p2^5sinp3 + PgP2, 




= Pi sin^ p3 - 2P3 COS p3 sin ps + /32 cos^ ps - 2pi P5 cos ^3 




+2;ji/34sinp3 + Pepj, 




= iPi- P2) sin p3 cos p3 + Ps (cos^ p3 - sin^ ps ) 




+(pi/34 + P2/35) cosp3 + {P1P5 - P2Pa) sinp3 - P&P1P2, 


Pa 


= Pi COS p3 + /35 sin P3 - PqP2 , 


P5 


= P^ COS p3 - /34 sin p3 - Pepi, 


Pe 


= Pe, 



(12) 



where and i = 1, ... ,6 are the arbitrary parameters that determine the factors in the 
product space P ([9]) according to the general formula (fTO]l . Recall that the group action 
SE{2) O /C^(E^) yields four types of non-trivial orbits (i.e., the orbits of dimesions 1, 2 and 3) 
corresponding to the elements of /C^(E^) with almost everywhere pointwise distinct eigenvalues. 
As is well known, these four types of orbits correspond to four types of Killing two-tensors whose 
geometric and algebraic properties are equivalent up to the action of the isometry group SE(2), 
namely the Killing tensors whose eigenvalues and eigenvectors generate elliptic-hyperbolic, po- 
lar, parabolic and cartesian coordinate systems respectively. Moreover, polar, parabolic and 
cartesian coordinate systems are all degeneracies of elliptic-hyperbolic coordinate system, char- 
acterized by the existence of two foci, which are the singlular points where the eigenvalues of 
the corresponding Killing two-tensor coincide (see, for example, |17j and the relevant references 
therein for more details). In our setting, namely for the classification problem determined by the 
group action SE{2) O -P, the analogue of the elliptic-hyperbolic coordinate system generated 
by a Killing two-tensor K £ /C^(E2), is a pair {Ki,K2) £ P, where Ki,K2 £ IC'^{E'^) are both 
of the elliptic-hyperbolic type. Moreover, the four foci of Ki and K2 are in the general position 
forming an arbitrary quadrilateral (with no symmetry structure). Let 81,82 and 53, 5*4 be the 
foci of Ki and K2 respectively. Denote by {xi,yi)si their respective coordinates given in terms 
of the cooresponding arbirtrary vector space parameters a's and /3's ()12p . Then, 



ix2,y2, 



S2 



^ \ 1/2 / ^ \ 1/2'' 

/36 /36 I 2 ) ' /36 /36 I 2 



/ ^ \ 1/2 / ^ \ 1/2^ 



(13) 
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53 



-"5 _|_ J_ 
"6 "6 



1/2 



-Q4 



"6 



^+0-2 



1/2N 



(14) 



Q6 



1/2 



-"4 
as 



1/2N 



where o"i = /3| — /3| + (^^{(^2 — Pi), (^2 = a\ — a1 + aQ{a2 — ai), while A3 and Ag are given below 

by (USD [iziia. 

According to the Fundamental Theorem on invariants of regular Lie group action [18] the 
action SE{2) O P given by (|12p yields 9 = dimP — dimS'£'(3) fundamental joint invariants. It 
was shown in [2] that the following 9 functions form a set of fundamental joint invariants of the 
group action (fT2]) : 





= tte, 








A^ 


= a6(ai + "2) 


- al ■ 








= {a&{ai - a2, 


) - al 


+ «i 


f + 4(a6a3 + oiAa^f 




= Pe, 










= P6{Pl+P2)- 


-Pi- 


-Pi 




K 


= {Pq{Pi-P2) 


-Pl 


+ Pl) 


.2+4(/36/33 + /34/35)', 


A'r 


= d^{S2,S3) = 


{X2 - 




+ (2/2 -2/3)^ 


A's 


= d^{Si,S3) = 


{xi - 




+ (2/1 -2/3)^ 


A'9 


= £{S2,S^) = 


{X2 - 




+ (2/2 -2/4)^ 



where d?{S2, S3), cf{Si, Ss), (f{S2, S^) are three joint ivariants representing the square distances 
between foci of Ki and K2. Consider the resulting quadrilateral S1S2SSS4. Obviously, it 
completely determines the geometry of the corresponding pair (i^i, 1^2) G P- Moreover, one can 
fix it by fixing the interpoint distances d{Fi, F2), d^S^, 84), d{S2, S3), d{Si, S3), d{S2, S/i), which 
are clearly invariants and joint invariants of the group action (I12p . Note that d{Si,S2) = 2ki, 
d{S3,Si) = 2k2, where 



2 _ V (al - "i + 06(02 - ai))^ + 4(a6Q;3 + a^a^)'^ 



,.2 _ ViP'i - Pi + PdP2 - Pi))^ + mPs 

Recall also that any joint invariant J {Si, 82,83, S4) of the non-transitive action so defined can 
be written as a function of the distances d{Si,S2), d{S2,S3), d{S3,S4,), d{Si,Si) (this result is 
also known as "the Weyl theorem on joint invariants" |18( I26j). 

Having constructed this groundwork, we can now proceed and use the language of the ITKT 
to invariantly characterize the Smorodinsky-Winternitz potential d?]). 
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4 Invariant characterization of the Smorodinsky-Winternitz po- 
tential 

As is well known the Smorodinsky-Winternitz potential of the Hamiltonian system defined by 
([7|) is superintegrable. Furthermore, its superintegrability is afforded by first integrals which are 
quadratic in the momenta ([4]), which follows immediately from the form of the SW potential. 
Indeed, since the potential in question is in the form V = f{x) + g{y) it follows that it is 
separable in (canonical) cartesian coordinates, which means that the Hamiltonian flow defined 
by ([Z]) admits a first integral Fi of the form: 



Fi = Kl'piPj + Ui, i,j = 1,2, (16) 

where = 1, Kf^ = Kj^ = Kf = 0. 

Next, since in the potential ([7]) is of the form V{r, 0) = + g{r) in polar coordinates, the 
Hamiltonian system defined by ([7]) is separable in (canonical) polar coordinates, which means 
that it admits the following quadratic first integral of motions: 

F2 = K'^'piPj + U2, i,j = 1,2, (17) 

where the components of the corresponding Killing two-tensor are as follows: iTg^ = y^, J^|^ = 
i^2^ = —xy, Er|^ = x^. It follows immediately that the Hamiltonian sysem defined by d?]) is 
superintegrable (note {Fi,F2} ^ 0) and multi-separable. It also follows that the SW potential 
is separable in (canonical) elliptic-hyperbolic coordinates, which is a consequence of the fact 
that the function F^ = Fi + F2 is a first integral of motion (see, for example, \n\ |2] for more 
details). Therefore, the Smorodinsky-Winternitz potential is separable in cartesian, polar and, 
- as a consequence, - in elliptic-hyperbolic coordinates. Separation in (canonical) polar and 
elliptic-hyperbolic coordinates implies that in terms of the corresponding joint invariants (see 
Section [3]) the quadrilateral S'i52<S'354 degenerates in this case. Let Si and S2 be the singular 
points (foci) of the Killing two-tensor of the first integral ^3 = ^1-1-^2) where Fi and F2 
are given by (fT6]) and pT|) respectively. 

Next, let S3 be the singular point of the Killing tensor defined by (I17p . Since, both K2 and 
K3 are in canonical form, we conclude that the quadrilateral S1S2S3S4 degenerates in this case, 
which is prompted by the invariant conditions S3 = S4, d(Si,S3) = d{S2,Ss), or, in terms of 
the fundamental joint invariants (jlSp . we have A3 = 0, Ag = Ag. 

Conversely, if we start with a general potential V given by ([T]) and then impose that the 
corresponding Hamiltonian system admits two quadratic first integrals of the form (|16p and 
(jl7p defined by Killing two-tensors Ki and K2 respectively, the corresponding compatibility 
conditions ([6]) given by Ki and K2 will yield the Smorodinsky-Winternitz potential ^ [H llOj. 

Conlcusion: the superintegrable system of the Smorodinsky-Winternitz potential is defined 
by a pair of Killing two-tensors {K2,K3) £ P, whose position in the orbit space P/SE{2) is 
determined by the following invariant conditions: 



a; / 0, A'3 = 0, A^ / 0, A'^ / 0, A'7 = A's = a;,. (18) 
Hence, we have 

Theorem 1. Let the potential V of the general Hamiltonian (OP be compatible via ^ with a 
pair of Killing two-tensors {Ki,K2) £ P- Then the following statements are equivalent: 

(1) The pair {Ki,K2) £ P is invariantly characterized by the conditions U8\) . 
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(2) The potential V given by lip is the Smorodinsky-Winternitz potential 

That is, by invariantly characterizing the pair (K2,K3) € P, we have thus characterized the 
Smorodinsky-Winternitz potential that they define. 

Now let us investigate what happens to the SW potential ([7]) if we weaken the invariant 
conditions (llSp . More precisely, let us first suppose that the Hamiltonian system defined by the 
natural Hamiltonian ([7]) admits two quadratic first integrals of motion 



Fi = Kl'piPj + Uu i,j = 1,2, (19) 

and 

F2 = Ki^PiPj + U2, i,j = 1,2, (20) 

Furthermore, we will assume that the Killing tensor Ki of (jl9p generates polar coordinates, 
while K2 of ()20p generates elliptic-hyperbolic and they are are in general position, which means 
their respective singular points ^3 (of Ki) and Si, S2 (of K2) form a general triangle AS'i5253. 
This arrangement entails the following invariant conditions in terms of the fundamental joint 
invariants (jlSp for the pair {Ki,K2) € P. 

a; / 0, A'3 = 0, A^ / 0, A^ / 0, A'j = A'g / A^. (21) 

Without loss of generality we can assume that K2 is in canonical form, while Ki is not, which 
means that the respective components of the two Killing tensors are as follows: 

Kl' = {y- h)\Kf = KP = -(x-a){y- h),Kf = {x-af,a,be R, (22) 

and 

= y^ + e, Kf = Kl^ = -xy, Kf = x^. A; G M, (23) 
where the parameter £ = 2k'^ is an invariant determined by the condition 

.2 _ \/(/54 - + W2 - /3i))2 + 4(/36/33 + _ 7(A^ 



/36 A^ 



Question: How does this setting affect the form of the Smorodinsky-Winternitz potential ([7|)? 



? 



Our first observation is that since -Fi and F2 are the first integrals of the Hamiltonian system 
defined by d?]), the Killing two-tensor given by the linear combination 

Kg = ciKi + C2i^2 + C3C/, ci, C2, C3 G M (24) 

must also be compatible with the SW potential V given by ([7]), i.e., d{KgdV) = 0. However, Kg 
depends on six arbitrary parameters, namely ^, a, b, ci, C2, C3, but recall that 6 is exactly the 
dimension of dim/C^(E^). This is only possible if Kg p4p is given by the general formula (llOp . 
which immediately implies that the potential ([7j) must be trivial, i.e., uj = a = j3. Therefore the 
only parameters that can be set to zero are a and b. Having discussed the case 

1) a/ 0,6/0, 

which reduced the SW potential ([7]) to a trivial one, now let us consider the following remaining 
three cases: 
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2) a = 0,6/0, 

3) a/ 0,6 = 0, 

4) a = 6 = 0. 

First, we note that in terms of the fundamental joint invariants these cases can be 
invariantly characterized as follows: 

1) a; / 0, A'3 = 0, A^ / 0, A^ / 0, A',A',A', / 0, A^ / A'g / A^, 

2) a; / 0, A^ = 0, A^ / 0, A'g / 0, the area of AS'iS'253 = 0. Note that in this case the 
singular points Si, S2, are on the x-axis. The area of AS1S2S3 is obviously also a joint 
invariant. 

3) A[ / 0, A3 = 0, A^i / 0, Ag / 0, A'j - A^ = 0. Note in this case AS1S2S3 is an isosceles 
triangle. 

4) a; / 0, A'3 = 0, A^ / 0, A^ / 0, A'7 = A'g / A^. 

Next, we note that the arbitrary parmeters a and 6 are obviously also joint invariants given 
in terms of the fundamental joint invariants ()15p by the following formulas: 



and 




(26) 



Solving the PDEs defined by the compatibility condition ([6]) for the Killing two-tensor (j22p and 
the SW potential ([7]), we arrive at the following equation 

2/L ^ , a{b-y)a b{-a + x)(3 
uj^{bx-ay)^ 2 \ = 0, (27) 

which puts in evidence that 

2) a = 0, 6 7^ yields the potential V = 

3) a 7^ 0, 6 = yields the potential V = 

4) a = 0, 6 = yields the general Smorodinsky-Winternitz potential ([7]). 

We conclude therefore that the more geometric structure the pair {Ki,K2) G P has, the 
less general the Smorodinsky-Winterniz potential ([7|) becomes. We also note that in each case 
described above the corresponding pair (7^1,^2) S P belongs to a different orbit of the group 
action SE{2) O P ()12p . Indeed, in the above the corresponding orbits have been distinquished 
by means of joint invariants (jlSp of the action (|12p . Extending in this paper the results obtained 
in [2], we have demonstrated that the invariant classification of the orbits P/SE{2) leads to a 
classification of the corresponding superintegrable potentials defined by quadratic first integrals 
in E2. 
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5 A remark on the Trembley-Turbiner-Winternitz potential 

Consider now the Trembley-Turbiner-Winternitz (TTW) potential defined by the natural Hamil- 
tonian ([8]). Clearly, it is separable in (canonical) polar coordinates, admitting a first integral of 
motion defined by the Killing two-tensor (I22p for a = b = 0. 

Question: For what value(s) of the parameter k does the TTW potential ([8]) admit an addi- 
tional functionally independent quadratic first integral of motion? 

Assuming k is real, we cannot transform the TTW system back into cartesian coordinates 
and so we must carry out all our calculations in terms of polar coordinates. Hence, we first 
rewrite the general formula the vector space (llOp in terms of polar coordinates as follows: 

Kg = A'^Xi Xj + B'Xi QR + l3eR, i,j = 1,2, (28) 

where is the symmetric tensor product and the generators Xi,X2, R of the Lie algebra se(2) 
are given by 



^ d sinO d , , 

^ d cos 9 d , , 

R = ^. (31) 



Here 

In view of the above the components of the general Killing tensors defining /C^(E^) in polar 
coordinates are given by 

K^^ = 2 133 cos e sin 9 + /3i cos 6^ + 132 sin 9"^, (33) 

12 2/33 COS 2g 2/3i cos g sing , 2/^2 cos g sing , . 
K = 1 h 2p4 cos 6 + 2/35 sm 0, (34) 



r^22 o 2/33 cos 6* sin g /3i sin 6*2 cos 9"^ 2/34 sin 6* 2/35 cos g 

K = Pq 5 1 5 1 5 1 . (35) 



Substituting (|28p and the potential of the natural Hamiltonian ([8|) into the compatibility 
condition ([6]), we arrive at a complicated trigonometric equation in terms of the parameters k, 
a, /3, uj, f3i, i = 1, . . . ,Q. To solve it, we expand the equation in terms of the following set M of 
trigonometric functions 

M = {l,cos(6l - 2mk),sin{9 - 2mk),cos{9 + 2mA;),sin(6' + 2mk),cos{29 - 2mk), 

(36) 

sin(26' - 2mA;),cos(26l 2mA:),sin(26' 2m/c)}, m = 0, ...,4. 

The set M is linearly dependent when the arguments of the trigonometric functions coincide, 
which depends on the values assumed by the parameter k. Indeed, we have that the set ([36]) is 
linearly dependent for the following values of k: 

k =±2,±l-,±l,±-,±-,±-,±-,± — ,± — ,± — ,±—,±-,±-, 

' 2' ' 2' 4' 6 8 10 12' 14' 16 4' 3' (oj^ 

3 1 3 2 3 1 3 2 1 ^ ' 

± -,±-,± — ,±-,±—,±-,±—,±-,±-. 
8 3 10 7 14' 5 16 5 7 
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In the case when k does not equal one of the above values, the set M is linearly independent and 
so we can set each of the coefficients to zero in the trigonometric equations that we obtained 
upon the substitution of (p8|) and ([8]) into ([6]). This argument leads to an overdetermined 
system of polynomial equations which we can solve for the variables k, a, /?, w, /3j, i = 1, . . . , 6. 
Though, the problem at this point can be simplified according to the following argument. Since 
the TTW potential ([8]) is separable in polar coordinates, if it were to admit another quadratic 
first integral of motion, it would have been defined by a Killing two-tensor which is an element 
of the 5-dimensional vector subspace 



i.e., the vector space (I28p with the "polar" term factored out. Substituting this Killing two- 
tensor Kg into the compatibility condition ^ for the SW potential we arrive at a system 
of equations, which is possible to expand over the basis L = {1, cos 9, sin 9, cos 29, sin 29}, which 
results in the following restriction on the parameters of the Killing two-tensor ()38p : 



Hence, we conclude that if the potential ([8]) is separable in a coordinate system, other than polar, 
it can only be the cartesian system of coordinates. Now we can repeat the same argument, but 
for a much simpler Killing two-tensor, namely the Killing two-tensor that defines cartesian 
coordinates in general position (not canonical) [l7], namely: 



where (j) denotes the angle by which the cartesian coordinate system is rotated. Substituting 
([39|) and ([8]) into the compatibility condition we get the following equation: 



Ks = A'^Xi Xj + B'Xi R, i,j = 1, 2, 



(38) 



Pi + /3i = 0. 




(39) 



- 4kp cos^{9 - (j)) cot{ke) csc'^ikO) + 4kp cot{k9) csc^{ke) sm^{9 - </>) 

- A^csc^{k9) sin(2(6' - 0)) + 2k^^cot^{ke) csc^{k9) sin(2(0 - 0)) 
+ k'^P csc^(A;0) sm{2{e - 0)) - 4/3 sec^{k9) sm{2{9 - (/))) 

+ k'^a sec^{k9) sm{2{9 - (p)) + 2k'^a sec'^{k9) sin(2(6' - (j))) tan^{k9) 
+ cos(2(6i - (p)) {-2k(3 cot{k9) csc'^{k9) + 6ka sec'^{k9) tan{k9)) = 0. 



(40) 



We can then expand this over a subset N C M given by 



N = {1, cos(2(l ± k)9), sin(2(l ± k)9)} 



m = 1, . . . , 3 



(41) 



which again is linearly dependent for the following values of k 



2 12 

k = ±l,±2,±-,±-,±- 
' ' 3 2' 5 



(42) 



which are indeed a subset of the values obtained in (j37p . Expanding the equation (j40p over the 
set N, and assuming it to be linearly independent (i.e. excluding the values of k given by (j42p ). 
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we arrive at the follwing system of polynomial equations in the variables (j), k, a, f3, namely 

{-1 + Ak^) + a)cos(f>sm(f) = 0, 

(-2 - k){l + /c)(/3 - a)cos</>sin(/) = 0, 

(1 + A:)(l + 2k){l3 + a)cos(j)sm(j) = 0, 

(-2 + k{l5 + 23/c))(/3 - a) cos (/)sin = 0, 

(-2 + k{-15 + 23A;))(/3 - q) cos sin = 0, 

(-1 + A;)(-l + 2/e)(/3 + a) cos sin = 0, 

(-2 + k){-l + /c)(/3 - a)cos0sin0 = 0, 

8 (-1 + 4A;2) (/3 + q) cos^ - 8 (-1 + ik^) (/3 + a) sin^ (j) = 0, 

(1 + A:)(2 + - a) cos^ + (-1 _ A;)(2 + - a) sin^ = 0, 

4/3 cos^ (/) + 12A:/3 cos^ (p + 8A;^/3 cos^ + 4a cos^ + 12A;a cos^ 4> 
+ 8A:^a cos^ (/> - 4/3 sin^ - 12A;/3 sin^ (p - Sk'^(3 sin^ cf) - Aa sin^ </> 

— 12A;a sin^ (j) — Sk'^a sin^ = 0, 

- 2/3 cos^ (/> - 15k(3 cos^ + 23k^(3 cos^ + 2a cos^ (/> + 15A;a cos^ (f> 

- 23A;^a cos^ (/) + 2/3 sin^ + l5k/3 sin^ - 23A;^/3 sin^ - 2a sin^ 

— 15A;a sin^ (j) + 23A;^a sin^ (j) = 0, 

4/3 cos^ (j) - I2kp cos^ + Sk^(3 cos^ + 4a cos (h — 12ka cos (/> 
+ 8A:^a cos^ (/> - 4/3 sin^ + 12A;/3 sin^ <^ - 8k^l3 sin^ - 4q sin^ (j) 
+ 12A;a sin^ — 8/c^a sin^ </> = 0, 

2/3 cos^ (p — 3kf3 cos^ (j) + /c^/3 cos^ (/> — 2a cos^ (j) + 3A;a cos^ 

- A;^a cos^ - 2/3 sin^ c/) + SkjS sin^ (j) - k^fi sin^ + 2a sin^ cf) 

— 3/ca sin^ cf) + A;^a sin^ = 0, 

(-2 + /c(15 + 23A:))(/3 - a) cos^ + (2 - A;(15 + 23A:))(/3 - a) sin^ = 0. 

There is no non-trivial solution for which the above system is satisfied. Considering the set of 
linearly dependent functions, we substitute into equation (j40p with the special values of k given 
in (jl2]) to yield the following possible cases 

1) k = ±1, 

= 0, which implies the polar coordinate system is in canonical form; 

2) k = ±2, ±2/3, ±1/2, ±2/5, 

w = a = /3 = 0, i.e., y is a trivial potential. 

Thus, we conclude that A; = ±1 are the only values of k for which the TTW potential dS]) admits 
two quadratic first integrals of motion. For these values of k it reduces to the SW potential ([7]). 
Hence, we have proven the following 



Proposition 5.1. The TTW potential given by ^ is a multi- separable superintegrale potential 
iffk = ±l. 
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6 Conclusions 

We have invariantly characterized the Smorodinsky-Winternitz superintegrable potential modulo 
the action of the isometry group SE{2) in terms of the joint invariants of the associated Killing 
two-tensors by treating them as points in the product space ([9]). This result demonstrates 
that it can be extended to the problem of classification of arbitrary superintegrable potentials 
defined in spaces of constant curvature. More specifically, such a classification can be carried 
out by classifying first any associated Killing tensors that define first integrals of motion of the 
class of superintegrable systems in question (e.g., the superintegrable systems that admit only 
quadratic in the momenta first integrals of motion). Such Killing tensors ought to be considered 
as elements of the corresponding product spaces, whose factors are the corresponding vector 
spaces of Killing tensors representing the corresponding orbits in the orbit space induced by the 
action of the isometry group in such product spaces. 

In addition, we have proven that the Tremblay-Turbiner-Winternitz potential is a superin- 
tegrable system with two functionally independent quadratic first integrals of motion only for 
/c = ±1, that is when it reduces to the Smorodinsky-Winternitz potential. 
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